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particle systems by exploiting diffusive processes 
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Over the past few years the displacement statistics of self-propelled particles has been intensely 
studied, revealing their long-time diffusive behavior. Here, we demonstrate that a concerted com¬ 
bination of boundary conditions and switching on and off the self-propelling drive can generate 
and afterwards arbitrarily often restore a non-stationary centered peak in their spatial distribution. 
This corresponds to a partial reversibility of their statistical behavior, in opposition to the above- 
mentioned long-time diffusive nature. Interestingly, it is a diffusive process that mediates and makes 
possible this procedure. It should be straightforward to verify our predictions in a real experimental 
system. 
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I. INTRODUCTION 

Diffusion .1] is one of the most studied and most fa¬ 
miliar phenomena in physics. Regular diffusion is well 
accessible on the different levels of description in statis¬ 
tical mechanics, making it a prominent textbook exam¬ 
ple On the single-particle level, the situation is 

typically modeled by corresponding Langevin equations 
00. These contain a stochastic force term to reproduce 
the experimentally observed Brownian motion. Likewise, 
the scenario can be characterized by a dynamic equation 
for the probability density to find a particle at a certain 
position at a certain time. This type of dynamic equa¬ 
tion is usually termed Smoluchowski equation or, more 
generally, Fokker-Planck equation [ 30 . The connections 
between the particle picture and these continuum equa¬ 
tions for the probability density are well established. Fi¬ 
nally, on a more macroscopic level, the phenomenon is 
addressed by monitoring the time evolution of concen¬ 
tration profiles using diffusion equations 130 . 

One particular feature of diffusive phenomena is their 
irreversibility. In general, it is not observed that an ini¬ 
tial density peak that has flattened over time will sud¬ 
denly rearise. This is reflected by the nomenclature in 
the theoretical approaches. One refers to the irreversible 
part of the Fokker-Planck operator 00 in the probabil¬ 
ity density picture and to the irreversible currents [3! in 
macroscopic descriptions. A deeper theoretical founda¬ 
tion is provided by the H-Theorem 00, which allows a 
connection to the increase of entropy or decrease in free 
energy. 

Here, we report on a globally diffusive system that is 
demonstrated to at least partially break this paradigm of 
irreversibility. Diffusive systems of initially peaked den¬ 
sity profiles generally do not reestablish such a density 
peak at a later time. Nevertheless, the example system 
analyzed below shows exactly such behavior, at least to 
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some degree. There, a non-stationary centered density 
peak can be generated independently of the initial distri¬ 
bution. Afterwards it can be restored arbitrarily often. 
It is the non-equilibrium nature of the investigated self- 
propelled particle system that allows for these effects, 
together with the possibility to switch on and off the 
self-propulsion. 


II. MODEL SYSTEM 

Self-propelled particles feature a mechanism of acti¬ 
vated motion [MU- On the one hand, this can for ex¬ 
ample be the crawling or swimming machinery of motile 
biological microorganisms such as bacteria and algal cells 
nang. On the other hand, self-propelled particles were 
artificially realized for instance in the form of granular 
hoppers fl9H22] or self-propelled droplets [2311251 . A par¬ 
ticularly interesting example are Janus particles that al¬ 
low localized heating by light illumination of only one 
half of their body [2H!l3()| : or they catalyze chemical re¬ 
actions on only part of their surface [3H - l33[ - The re¬ 
sulting build-up of temperature or concentration gradi¬ 
ents along their surfaces induces stress differences that 
can effectively drive the particles forward in a type of 
phoretic motion [[34] . As a central common property, the 
migration direction of self-propelled particles is not im¬ 
posed or fixed from outside but results from symmetry 
breaking [10]. Light-activated self-propulsion offers the 
advantage of being easily switchable by turning on and 
off the external illumination PBM3U1 [551 . 

We consider for our purpose mostly the statistics of the 
motion of non-interacting self-propelled particles. This 
implies that either dilute systems are investigated, or the 
averaged statistical behavior of single objects is analyzed. 
At the end we will briefly address the effect of basic steric 
interactions. As is natural e.g. for microswimmers, we ad¬ 
dress their overdamped motion. For simplicity and since 
many experiments are performed in this way [I5l - l24| l26t - 
132] , we confine ourselves to two spatial dimensions. The 
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corresponding Langevin equations for a single particle 
then read (compare e.g. Refs. 32^ 1361) 

dr 

- = (3D t [F u - W(r)] + y/2 A €„ (1) 

= \j2D r ^. ( 2 ) 

Here, r denotes the particle position and u the orienta¬ 
tion of its self-propulsion direction. F sets the strength 
of self-propulsion that drives the particle forward along 
u. The potential V(r) leads to a force acting only on 
the particle position. In our case, it will be nonzero only 
at the system boundaries to include the steric effect of 
a confining cavity, t denotes time and /? = (fcsT )^ 1 
with ks the Boltzmann constant and T the tempera¬ 
ture. In the two-dimensional plane, the self-propulsion 
direction can be parameterized by a single angle 0 via 
u = (cos 0, sin 0). D t sets the translational and D r 
the rotational diffusion constant. Finally, £ r and 
are Gaussian white noise terms with averages (£ r (t)) = 
0 , (&(£)> = 0, <&,<(*)&(*')> = 0, <&,<(*)&,i(0> = 
6ijS(t — t'), and = 5(t — t'), where Sij is the 

Kronecker delta and 5(t — t') the Dirac delta function. 

We now switch to the probability density description 
by deriving from Eqs. ([T]) and ([ 2 ]) the corresponding 
Smoluchowski or Fokker-Planck equation [2 0133. In 
dimensionless units, this dynamic equation for the prob¬ 
ability density 0(r, £) to find a particle at time t and 

position r with orientation 0 of the self-propulsion direc¬ 
tion becomes 

dttp = — vq (cos 0 d x ip + sin 0 d v ip) + V 2 0 + 3^0 

+ V • [(W)0]. (3) 

Here, we rescaled time, space, and the potential by 
t! = tD r , r' = YyjDr/Dt, and V' = /3V, respec¬ 
tively, with the primes omitted in the equation, and we 
confine ourselves to the x-y plane. The only remain¬ 
ing parameter determining the system behavior is now 
vo = /3F\jD t lD r , which measures the relative strength 
of the self-propelling drive. From the data listed in 
Ref. [57] we infer that vq = 30 constitutes a reasonable 
order of magnitude. 

In our case, the potential V (r) only contains the effect 
of the cavity boundaries. Within the bulk of the cav¬ 
ity, the potential vanishes and the particle is completely 
free. Furthermore, V does not act on the orientation of 
the self-propulsion direction u. This direction is solely 
determined by the independent rotational diffusion pro¬ 
cess Eq. <§• It has been demonstrated in various works 
that, as a consequence, the global and long-time transla¬ 
tional behavior of such particles is itself diffusive, with an 
effective diffusion constant increased by self-propulsion 
[50 [571 03 [30. Despite these facts, we identify in the 
following a procedure that breaks the irreversibility im¬ 
plied by this long-time diffusive behavior, at least par¬ 
tially. Remarkably, our procedure is mediated by and 
only becomes possible due to rotational diffusion as de¬ 
scribed by Eq. ([2]). 


III. CHANNEL 

We start by considering an infinitely extended chan¬ 
nel. This implies two parallel and infinitely extended 
confining walls, in our case along the y direction. For 
simplicity, we choose the confining potential as 

V'(r) = V(x) = { k (1*1 “ ^ lf M>T’ ( 4 ) 

1 0 if |*| < f, 

where the width of the channel between the wall regions 
is given by L > 0. The precise magnitude of the exponent 
and of the potential strength k > 0 is not qualitatively 
important for the presented results. We set k = 50. 

As an initial condition, we consider self-propelled 
particles located in the center of the channel with 
equally distributed orientations u. The correspond¬ 
ing probability density is independent of the y coordi¬ 
nate, 0(a;,0, f = 0) = S(x)/2tt. Consequently, the prob¬ 
lem becomes translationally invariant in the y direction. 
The spatial probability profile 'tj) x (x, t) across the chan¬ 
nel is obtained by integrating out the 0 -dependency of 
i/j(x,(f>, t). We choose a channel width of L = 20, which 
in the framework of Ref. m corresponds to about 11 
particle diameters. 

Our protocol is the following. We numerically integrate 
the resulting dynamic equation for i/j(x, 0 , t) forward in 
time. Fig. [l] (a) shows the sharply-peaked initial proba¬ 
bility density profile il> x (x,t ~ 0 ) at a very early stage, 
together with the confining potential V{x). Due to self¬ 
convection and since all orientations of the self-propulsion 
direction are equally probable at t = 0 , this peak splits 
and two fronts propagate outward towards the channel 
boundaries, see Fig. |T] (b). At the boundaries, the prob¬ 
ability density piles up and after a while reaches a sta¬ 
tionary profile as displayed in Fig. [l] (c). Such a blocking 
effect due to a confining potential has been observed in 
various forms before [571 EEH10 : the self-propelling drive 
pushes the particles towards the walls, and the only way 
of escape is a reorientation of the self-propulsion direc¬ 
tion. In our case, this reorientation is solely due to rota¬ 
tional diffusion. 

We now switch off self-propulsion, setting vq = 0. As 
mentioned above, for light-activated mechanisms this can 
simply be achieved by turning off the external illumina¬ 
tion (3003 S0. This has two effects as displayed in 
Fig. [l] (d). First, the potential pushes the particles back 
towards the edges x = ±L/2. And second, translational 
diffusion leads to a broadening of the density peaks. Nev¬ 
ertheless, rotational diffusion is still active and generates 
within this localized particle cloud equally strong sub¬ 
populations of inward and outward oriented particles. 

When then self-propulsion is switched on again, v 0 7 ^ 0, 
the peaks at the boundaries split. Outward oriented par¬ 
ticles return towards the channel boundaries. However, 
inward oriented particles form peaks that start to propel 
back towards the channel center, see Fig. l] (e). When 
they meet in the channel center, they overlay and again 
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FIG. 1. Time evolution of the spatial probability density 
distribution i/j x (x,t) across the channel of width L, bounded 
by a confining potential V ( x ) of strength k. From the initial 
density peak (a) fronts migrate outward (b) and get blocked at 
the confining boundaries (c), where stationary density peaks 
build up. Switching off the self-propelling drive, these peaks 
broaden due to spatial diffusion (d); reorientation is possible 
due to rotational diffusion. Turning on the drive again (e), the 
peaks at the boundaries rebuild, but also two further peaks 
start to head back towards the channel center. When they 
overlay in the channel center (f), the initial spatial density 
distribution is partially restored. From here, the cycle can be 
repeated arbitrarily often without further losses. Dominating 
self-propulsion directions are marked by arrows. Parameter 
values are L = 20, Vo = 30, and k = 50 in rescaled units. 


FIG. 2. Detailed probability density distribution tp(x, <p, t ) 
corresponding to Fig. [T| (b)-(e) in the x-(p plane. At early 
times (a), the initial front at x = 0 has split, leading to a 
peak at <p = 0, 27r that travels to the right, while the one 
at (j> = 7r migrates to the left. The peaks get blocked and 
build up at the confining boundaries (b), centered around 
the self-propulsion directions that point outwards. When the 
self-propelling drive is switched off (c), these peaks broaden, 
leading to fronts of equally distributed orientations (/>. Fi¬ 
nally, when self-propulsion is turned on again (d), those parts 
of the fronts corresponding to outward orientations reform 
the blocked peaks of (b); those parts of the fronts containing 
inward orientations start to propagate towards the channel 
center. Arrows mark the inward and outward self-propulsion 
directions. 


form a (smaller) central peak, as shown in Fig. 0(f)- 
Consequently, the initial distribution has been restored 
partially, with significantly lower magnitude of the cen¬ 
tral peak (we show below that the restored peak height 
can be notably increased by steric interactions). Af¬ 
terwards, following the same protocol, this central peak 
can be restored arbitrarily often without further losses. 
The restoration is possible despite the irreversible diffu¬ 
sive processes involved. In fact, we have explicitly ex¬ 
ploited a diffusive process to achieve our goal, namely 
rotational diffusion during the period of switched-off self¬ 
propulsion. 

Fig.0shows the whole process in the x-<f> plane. There, 
we explicitly observe the initial splitting into outward 
propagating peaks centered around <f> = 0, 27t and <j> = 7r, 
see Fig. [2] (a). The outward oriented particles get trapped 
at the boundaries, leading to the peaks localized both in 
x and </> direction in Fig. [ 2 ] (b). After switching off the 
drive, rotational diffusion leads to equally ^-distributed 
localized states at the boundaries. In other words, the 
peaks broaden to walls extended in the </> direction in 
Fig. i (c). Finally, when the drive is turned on again, 
the inward oriented parts of these walls start to propagate 
back to the channel center, while the outward oriented 
ones return towards the wall regions, see Fig. 0(d)- 

During this cycle, it is necessary to switch off the 


drive to generate with sufficient probability particle ori¬ 
entations towards the channel center. Turning on the 
drive, such particles immediately head towards the cen¬ 
ter, forming the propagating fronts that rebuild the cen¬ 
tral peak. If the drive were not switched off, a parti¬ 
cle at a wall would always leave the boundary as soon 
as it slightly turns away. At this moment, due to the 
typically gradual rotational diffusion, the main com¬ 
ponent of the drive points along the bounding surface 
and not towards the channel center; pronounced inward- 
propagating fronts could not form. 

The time scales for the individual processes are de¬ 
termined by the intrinsic system parameters. It takes 
the ballistic time L/2vq ~ 0.33 for the edge of the out¬ 
ward propagating fronts in Fig. 0 (b) to reach the chan¬ 
nel boundaries. One should wait an order of magnitude 
longer for the stationary peaks in Fig. 0 (c) to build up. 
Then, after switching off the drive, rotational diffusion 
over a time span At leads to a mean squared angular dis¬ 
placement of ((A0) 2 ) = 2Af in rescaled units. To achieve 
7PF) > 7T to smoothen out the orientational distri¬ 
bution, a period of At > 7 t 2 /2 ss 5 is necessary. Yet, 
one should not wait too long so that the peaks are not 
smeared out by translational diffusion, see Fig. 0(d)- 
Roughly estimating the mean squared displacement for 
the one-sided diffusion by ((Ax) 2 ) ~ At in rescaled units 
and requiring yj ((Ax) 2 ) < Lj8 , we obtain an upper con- 
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FIG. 3. Same process as in Fig. [I] but now without the 
confining channel walls, i.e. V(x) = 0. Spatial distances are 
rescaled by the same length L for comparison. From the ini¬ 
tial density peak (a) fronts migrate outward (b). Switching 
off the self-propelling drive, these free-standing peaks broaden 
due to spatial diffusion (c); reorientation is possible due to ro¬ 
tational diffusion. Turning on the drive again (d), the peaks 
split into two counterpropagating fronts. When the inward¬ 
traveling peaks overlay in the center (e), the initial spatial 
density distribution is partially restored. Over time, the den¬ 
sity profile continuously flattens (f) due to the absence of 
confining boundaries. 


straint for the waiting time At < L 2 /64 ss 6. After turn¬ 
ing on the drive again, the ballistic time L/2vq ss 0.33 
determines when the peak in the center of the channel 
reforms in Fig. |T] (f). To observe the formation of the 
central peak, the channel center must be within reach of 
the self-propelling front before it is diluted by diffusion, 
i.e. L /2 < vq. 

We should mention that a similar process could also 
be achieved without boundaries. In Fig. [3] we switch off 
the drive after the initial separation, let the rotational 
diffusion work, and then switch on the drive again. This 
likewise leads to an overlay and partially restored peak 
at the initial starting point. However, over time, the 
peaks will more and more broaden due to diffusive pro¬ 
cesses, and density is lost towards the open boundaries. 
In contrast to that, a stationary state is involved in the 
form of the built-up peaks at the boundaries in Fig.[l](c). 
From here, the central peak in Fig. [l] (f) can in principle 
be generated in an infinite number of repeatable cycles 
without further loss of restored peak intensity. The pro¬ 
cess becomes independent of initial conditions. Intervals 
of self-propulsion vq ^ 0 are necessary for this process to 
be observed. 


IV. CIRCULAR GEOMETRY 

In an experiment using a channel structure, the chan¬ 
nel ends may influence the statistics. Therefore, we addi- 


FIG. 4. Time evolution of the spatial probability density dis¬ 
tribution ip r (r, t) in the circular cavity of diameter L, bounded 
by a confining potential V(r) of strength k. From the initial 
central density peak at r = 0 (a) a front migrates outward (b). 
It gets blocked at the cavity boundary (c), where a station¬ 
ary density peak builds up. When the self-propelling drive is 
switched off, this peak broadens due to spatial diffusion (d); 
simultaneously, the orientations become equally distributed 
within the peak due to rotational diffusion. Turning on the 
drive again (e), the peak at the boundary rerises; an addi¬ 
tional second peak emerges and heads back towards the cavity 
center. When it arrives there (f), the initial spatial density 
distribution is partially restored. From here, the cycle can be 
repeated arbitrarily often without further losses. Dominating 
self-propulsion directions are marked by arrows. Parameter 
values are L = 20, Vo = 30, and k = 50 in rescaled units. 


tionally turn to a finite confinement. We consider a two- 
dimensional circular cavity of the same diameter L as 
the channel width above. The confining potential V (r), 
with r = Ill'll, has the same form as V(x) in Eq. 0- 
We assume an initial density distribution peaked in the 
center of the cavity at r = 0, again with equally dis¬ 
tributed orientations of the self-propulsion direction u. 
Thus the whole process becomes radially symmetric. We 
can reduce the number of variables measuring the self¬ 
propulsion direction by an angle i9 relatively to the radial 
direction. For the resulting probability density ij}(r,'d,t) 
we obtain 


d t 1p = — Vo COS 1? drip -sin 1? dfilp 

r 


- dr [r d r ip] + 

r 


1 + ^2 


d#tp + -d r [r(d r V)ip ]. (5) 


We now follow the same protocol as before. The pro¬ 
cess is displayed in Fig. [4] for the radial probability den¬ 
sity ipr( r it) obtained by integrating out the angle d in 
ip(r, d, t). Starting from a density peak around r = 0 
in Fig. [4] (a), particles migrating outward from the cen¬ 
ter of the cavity as in Fig. [I] (b) get trapped at the 
boundary. There, a localized peak of elevated proba¬ 
bility density forms, see Fig. [4] (c). Switching off the 
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FIG. 5. Same system as in Fig. [l] for the channel geometry, 
but now with an additional steric pair interaction v(r, r') = 
eS(r-r'). For increasing interaction strength e, starting from 
t = 0 , (a) the height of the stationary boundary peaks for 
Vo 7 ^ 0 decreases, whereas (b) the height of the restored central 
peak increases. The peaks are smeared out at elevated values 
of e. 


FIG. 6 . Same system as in Fig. [4] for the circular geometry, 
but now with a steric pair interaction u(r, r') = eS(r — r') 
included. Starting from e = 0, (a) the height of the stationary 
boundary peaks for vo ^ 0 decreases, while (b) the height of 
the restored central peak increases for increasing interaction 
strength e. At pronounced magnitudes of e the peaks are 
smeared out. 


self-propelling drive, rotational and spatial diffusion re¬ 
main active and lead to equally distributed orientations 
within the broadening density peak, see Fig. [4] (d). When 
the self-propelling drive is turned on again in Fig. [4] (e) 
inward-pointing particles start to head towards the cav¬ 
ity center. After a while, a density peak is observed in 
the center of the confinement in Fig. [4] (f). As before, the 
initial distribution has been partially restored. Further¬ 
more, the situation in Fig. [4] (f) can afterwards be com¬ 
pletely restored arbitrarily often and does not depend on 
the initial condition any more. This is because an inter¬ 
mediate stationary state was involved in Fig.[4](c). 

V. STERIC INTERACTIONS 

In Refs. [2HH2E], statistical properties of isolated light- 
switchable self-propelled particles were experimentally 
determined for dilute systems. We now briefly discuss 
what happens to the above results when particle interac¬ 
tions start to become important. 

The nature of steric mi m as well as of possible 
hydrodynamic interactions [431 ill] depends on details 
of the particular system under consideration. To re¬ 
main at this stage as general as possible, we only con¬ 
sider the most central feature of hard steric pair interac¬ 
tions: a diverging interaction energy when two particles 
are found at the same position. The simplest way to 
include this general property is by a steric pair inter¬ 
action of the shape v(r,r') = ei5(r — r'), where r and 
r' mark the positions of the two particles and e is the 
rescaled interaction strength. In the derivation of the 
statistical equations, we apply the mean-field approx¬ 
imation for the two-particle density by inserting 
■i// 2 )(r, cj), r', ft, t) = i />(r, 0, t)i/>(r', <(/, t) [43]. This leads 
to an additional contribution eV- [^(r, <^>)V /'0(r, cf)')d(j)'] 
in Eq. 0. 

There is now a certain pressure on the particles to leave 
highly populated areas. It is reflected by a lower mag¬ 
nitude of the stationary boundary peaks for vq 0, see 
Figs. [5] (a) and [6] (a) for the two considered geometries, 
respectively. Remarkably, close to the dilute limit, where 


the above approximations are restricted to, the steric in¬ 
teractions as a result enhance the build-up of the restored 
central peak, see Figs. [5] (b) and[b] (b). At higher inter¬ 
action strengths, the density distributions are smeared 
out. 

VI. CONCLUSIONS 

In this work, we have demonstrated how in a confined 
self-propelled particle system a central density peak can 
be repeatedly generated. Our procedure implies switch¬ 
ing on and off the self-propulsion mechanism in combina¬ 
tion with appropriate confining boundary interactions. If 
we start from an accordingly generated (or otherwise pre¬ 
pared) initial central density peak, this initial peak can 
(partially) be restored. After the first cycle of restora¬ 
tion, the central density peak can be repeatedly restored 
arbitrarily often without further losses in peak magni¬ 
tude. In that sense, we break the irreversibility of this 
long-time diffusive system. 

Naturally, there are other examples where diffusive 
processes can be reversed. For instance, the diffusion 
of dipolar electric or magnetic particles can be counter¬ 
acted by a guiding drift force resulting from an electric 
or magnetic field gradient. Yet, in such situations, the 
external forces directly pull the particles, imposing and 
pre-setting the drift directions. The present case is dif¬ 
ferent. Here, each particle itself individually selects its 
migration direction in a stochastic process of rotational 
diffusion. 

Many experiments on self-propelling particles are per¬ 
formed in two-dimensional set-ups P3H2H G6HH3. Ac¬ 
cordingly, we here concentrated on two spatial dimen¬ 
sions. In principle, the analogous effect could also be 
analyzed in three dimensions. However, this renders the 
analysis considerably more tedious and requires a well- 
resolved three-dimensional detection technique on the ex¬ 
perimental side. 

It should be straightforward to verify our results in 
corresponding experiments. Switching on and off the 
self-propulsion can easily be achieved for light-activated 
mechanisms [2BH50I [55 ]. As briefly demonstrated, steric 
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interactions among the particles can promote the ef¬ 
fect. Details of interactions between particles maim edi, 
steric alignment interactions between non-spherical par¬ 
ticles and cavity walls [TUI [HI GUI or hydrody¬ 

namic interactions among the particles and with cavity 
walls HOI HU HU HU HS1 HU HU may influence the statis¬ 
tics in different ways. The above theoretical description 
should then be adjusted to the specific situation under 


investigation. 
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